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Abstract 

We determine the bottom quark mass from non-relativistic large-n T sum rules with renormaliza- 
tion group improvement at next-to-next-to-leading logarithmic order. We compute the theoretical 
moments within the vNRQCD formalism and account for the summation of powers of the Coulomb 
singularities as well as of logarithmic terms proportional to powers of a s ln(n). The renormalization 
group improvement leads to a substantial stabilization of the theoretical moments compared to pre- 
vious fixed-order analyses, which did not account for the systematic treatment of the logarithmic 
a s ln(n) terms, and allows for reliable single moment fits. For the current world average of the strong 
coupling (a s (M z ) = 0.1183 ± 0.0010) we obtain M^ s = 4.755 ± 0.057 port ± 0.009 Qs ± 0.003 cxp GeV 
for the bottom IS mass and m^m^) = 4.235 ± 0.055 pert ± 0.003 exp GeV for the bottom MS mass, 
where we have quoted the perturbative error and the uncertainties from the strong coupling and 
the experimental data. 
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I. INTRODUCTION 



Determinations of the bottom quark mass m b have a long history and belong to the 
early precision analyses carried out in QCD that required higher order computations in 
the framework of the operator product expansion (OPE) [T] [2]. Phenomenologically they 
are an important input in the theoretical description of semi-inclusive B-meson decays and 
studies of T-mesons, but they also have some relevance for Higgs decays into 6-jets or 
unification and beyond the standard model studies. Apart from their phenomenological 
relevance, precision determinations of heavy quark masses are also a reflection of the ability 
to do precision computations for the quantities from which they are extracted, so that 
continued work in refining the respective theoretical methods and cross checks with other 
determination methods are warranted. Since quark masses are not physical observables due 
to the confinement, which binds quarks into hadronic states, they represent renormalization- 
dependent formal field theoretic quantities that are defined through the renormalization 
conditions imposed to remove the ultraviolet divergences in loop computations. Since the 
achievable precision in the bottom quark mass is at the percent level, and thus substantially 
smaller than the hadronization scale Aqcd ~ 200 MeV, it is mandatory to adopt so-called 
short-distance mass schemes for the bottom quark mass, which are free of O (Aqcd) infrared 
renormalon ambiguities. Short-distance mass schemes involve infrared subtractions in order 
to reduce the sensitivity of the perturbative series to small momenta. Thus short-distance 
masses either have an intrinsic or an explicit dependence on the infrared subtraction scale 
which should be compatible with the typical short-distance scale governing the heavy quark 
mass-dependent observable [3H3]. 

One of the classic methods to determine the bottom quark mass are the T sum rules. They 
involve moments P n of the total bottom pair production cross section in e + e~ annihilation 
through a virtual photon, defined by 



where R bb = a(e + e~ — > bb + X)/a pt is the inclusive bottom pair production cross section and 
Opt = 47ra 2 /3s is the point cross section. Conceptually the moments P n are related to the 
coefficients of the zero-momentum expansion of the photon vacuum polarization function 
coming from bh quark pairs [B] 



which can be determined from an OPE where the power series is dominated by the leading 
perturbative contribution and the power corrections represent an expansion in Aqcd over 
the short-distance scale m b . 

Since for large values of n the energy range AE, where R bb can effectively contribute to P n 
in Eq. (Til) , scales as m b /n [7], the scale m b /n emerges as an additional short-distance scale 
in the OPE and the values of n one can use for reliable perturbative analyses are limited up 
to around 10. Concerning the methods required for theoretical computations, two regions in 
n can be distinguished. For n < 3 the theoretical moments are dominated by fluctuations at 
the scale m b , such that the usual loop expansion is sufficient to determine the perturbative 
contributions. A suitable short-distance mass for the low-n T sum rules is the MS mass which 
has the bottom mass itself as the intrinsic infrared subtraction scale. On the experimental 
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side, the continuum region for R h i has a substantial contribution to the moments. Since 
there are sizable energy regions above the T bound state regime where there aren't any 
experimental measurements, one has to rely on extrapolations for the experimental moments 
P^ xp - Analyses available in the literature for the low-n sum rules differ concerning the 
uncertainty associated to this extrapolation, see e.g. the analyses in Refs.jS] and [9J. There 
are variants of this method, called finite-energy sum rules where certain energy ranges are 
excluded entirely from the analysis [ID] ■ For finite-energy sum rules systematic methods for 
computing power corrections in the OPE are currently unknown. 

For 4 < n < 10 the moments are dominated by the T bound state region and non- 
relativistic heavy quark pair bound state dynamics. The large-n sum rules are therefore 
also called non-relativistic or T sum rules. They are the focus of the present work. For the 
large-n sum rules, in addition to rrib, also m^j yfn and m^jn emerge as additional scales in 
association to the hard, soft and ultrasoft momentum scales known from heavy quarkonium 
physics, and one has to rely on non-relativistic effective field theory methods based on 
NRQCD [UJ to determine the perturbative contributions. Suitable mass definitions belong 
to the class of the so-called threshold masses, which have the soft scale rrib/y/n as the intrinsic 
infrared subtraction scale. Several threshold masses have been defined in the literature [2H5l 
IT2] . Concerning experimental data, no extrapolations are required, since the experimentally 
very well measured T resonances entirely dominate the moments P^ xp and the continuum 
contributions of R b i where no measurements exist are strongly suppressed. For a summary 
of various bottom mass determinations we refer to Ref. [TB"] . 

In this work we are concerned with the large-n non-relativistic sum rules. Because of the 
additional scales, perturbation theory is a simultaneous expansion in a s and 1/ '\/n, where 
terms proportional to powers of (a s \/n) or (a s ln(n)) count as order unity and are summed 
to all orders. The Coulomb singular terms oc (a s \/n) k are generated from ladder-type 
diagrams with iterated potential interactions and in the so-called fixed-order expansion only 
these terms are being summed systematically within the NRQCD formalism [llj. Large-n T 
sum rules analyses were first considered in the fixed-order expansion by Voloshin in Ref. [7] at 
leading-order (LO) and improved up to next-to-next-to-leading order (NNLO) in Refs. [T4T 
ITS] with differing treatments of the logarithmic terms and the residual higher order terms and 
with different choices of threshold masses. Unfortunately, the moments in the fixed-order 
expansion do not exhibit any good convergence properties and in Refs. [X4TTI7] various ad hoc 
methods were devised to extract the bottom quark mass at NNLO with a small uncertainty. 
We refer to the review [19] for a comparison and comments. Although the results of these 
fixed-order analyses were compatible, the situation was certainly unsatisfactory. We will 
comment on results obtained in the fixed-order expansion in Sec. |Vj 

In the renormalization group improved (RGI) expansion also the logarithmic terms are 
summed systematically with the order counting scheme 

P n ~ ^(a sV ^) fc KlnH)' [l(LL); a s , 1/ \/n (NLL); a 2 s , a s / \/ri, 1/n (NNLL); ...] (3) 

k,l 

for the leading logarithmic (LL), next-to leading logarithmic (NLL) and next-to-next-to 
leading logarithmic (NNLL) order. The RGI computations can be carried out in extensions 
of the original NRQCD formalism. The pNRQCD formalism [20J is based on the static 
expansion and factorizes soft and ultrasoft modes. The vNRQCD formalism [21j is based on 
the correlation of soft and ultrasoft modes through the energy-momentum dispersion relation 
for dynamical heavy quarks. Both formalisms are orthogonal concerning their conceptual 
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basis and have different Lagrangian formulations, but results which became available within 
both formalisms agree. 

Prior to this work a RGI bottom quark mass determination in the RS-scheme was carried 
out in Ref. [22J (see also [23J) at NLL order. This analysis also contained NNLL order 
matrix element corrections, but not the NNLL order corrections to the renormalization 
group evolution of the dominant quark pair production current which directly affect the 
normalization of the moments P n . The NNLL order corrections to the current anomalous 
dimension have a contribution from genuine three-loop diagrams, called non-mixing terms, 
and a contribution from the subleading running of the Wilson coefficients entering the NLL 
anomalous dimension, called mixing terms. The ultrasoft contributions to the subleading 
running of the relevant coefficients were determined recently in Refs. [24"1 [25] in vNRQCD 
and in Ref. [26J in pNRQCD. Although the computations differed substantially concerning 
their conceptual setup and computational effort, the results from vNRQCD and pNRQCD 
agree. In Refs. [MJ |25] also the complete set of the NNLL ultrasoft mixing contributions 
to the current anomalous dimension were given. The complete set of non-mixing terms has 
been determined already some time ago in Ref. [27] in vNRQCD, and for the consistent 
computation the correlated treatment and the simultaneous presence of soft and ultrasoft 
modes was essential. This is due to a term oc hi(/z| / 'rribHu) appearing in the anomalous 
dimension, which enforces the relation fin oc fi^/mj, between the soft renormalization scale 
Us and the ultrasoft renormalization scale The corresponding calculations have not yet 
been achieved in pNRQCD. The results for the mixing and non-mixing terms each exhibit 
very large corrections from ultrasoft gluons. As was found in Ref. [23] , these large ultrasoft 
contributions cancel to a large extent in the sum, confirming a speculation made in Ref. [27J. 
Also, the known soft corrections to the current anomalous dimension, which include the 
complete set from the non-mixing terms [27] and the spin-dependent contributions from 
the mixing terms [28], are much smaller than the ultrasoft ones and even smaller than the 
typical residual scale dependence of the combined ultrasoft corrections [29J . 

In this work we determine the IS bottom quark mass [5] from the large-n moments 
P n in the RGI expansion including all known NNLL order corrections, in particular the 
new contributions to the Wilson coefficient of the dominant quark pair production current. 
Our analysis is based on fits of single moments, and we examine in detail the convergence 
properties of the moments and their perturbative uncertainties. For the still unknown NNLL 
soft mixing contributions to the anomalous dimension of the dominant quark pair production 
current we give arguments that their size is comparable to the known soft contributions. The 
resulting uncertainty is tiny, and in particular smaller than the remaining scale variation 
of the known ultrasoft corrections. We therefore consider our analysis as NNLL order. 
Compared to the situation in the fixed-order expansion we find a reasonably good behavior 
of the perturbation series which does not require any ad hoc methods to determine the 
bottom quark mass. 

The outline of this paper is as follows: In Sec. [IT] we discuss the theoretical input for the 
non-relativistic factorization formula for the NNLL order moments we use for our analysis. 
For the most part we review older results, but we also discuss and numerically analyze 
the new NNLL order corrections to the anomalous dimension of the leading quark pair 
production current. An important aspect of the numerical analysis is that these NNLL order 
corrections are sizable. In Sec. Ill we therefore analyze the perturbative expansion and find 
that only fully expanding out all terms in the NNLL order factorization formula leads to the 
desired non-relativistic renormalization scaling behavior. Our numerical analysis concerning 
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fits for the bottom mass determination and the details on how we implement renormalization 
scale variations and estimate the perturbative uncertainty are discussed in Sec. IV In Sec. [V] 
we discuss the results of the fits when we modify our moments to account only for the 
terms required in the fixed-order expansion and demonstrate the improvement achieved 
by the resummation of logarithms. We also present arguments that disfavor carrying out 
simultaneous fits with several moments using a \ 2 function due to the strong correlation 
between different large-n moments. Section |VI| contains the conclusions. We have also 
added three appendices giving details on the analytic results for NNLL order moments and 
presenting the experimental moments. 



II. THEORY MOMENTS WITH RG IMPROVEMENT 

For the determination of the theoretical expressions for the large-n moments with renor- 
malization group improvement we follow closely the approach of Refs. 0[TH[T7] to implement 
the simultaneous expansion in a s and l/\/n according to the scheme in Eq. and, here, 
mainly concentrate on explaining the outline of the calculations and the new aspects that 
arise for the RG-improved calculation. The NNLL order large-n moments can be determined 
from the integral 

oo 

pth.NNLL = 1 f dE f_ E n \(l_ E i E<1 n \ R thv (E) 

'^{Ml ole ) 2n J Ml ole I Ml ole A 2Mf ie 4(Mf' e ) 2 / N 

(4) 



£1 



bind 



where E = <Js - 1M v h ole is the energy w.r.t. the pole mass threshold and -Ebind is the 
binding energy of the lowest lying IS resonance. Compared to the original definition of 
Eq. ([!]) the weight function l/s n+1 has been rewritten as an exponential expansion that 
explicitly separates the LL contributions from the NNLL order kinematic corrections. For the 
normalized bb total cross section i?-ratio in the threshold region at NNLL order, R% hh , one 
has to insert the absorptive part of the leading order and 0(v 2 )-suppressed non-relativistic 
vector current correlators arising for the heavy quark pair production process through a 
virtual photon in e + e~ annihilation. The required results are well known from the literature 
on top pair production in the threshold region at a future linear collider and can be applied to 
the bottom pair cross section with trivial adaptations mostly related to the different number 
of active flavors. The explicit results can be taken from Ref. [30H32] supplemented by the 
modifications given in Ref. [27] arising from a more convenient convention for the l/(mk)- 
type potentials. The contributions to the moments related to the Coulomb interaction 
are the same as for the fixed-order moments up to the additional logarithms summed in the 
leading order Coulomb potential Wilson coefficient and have been given explicitly in Ref. [H] . 
Deforming the energy integration into the complex plane the integral in Eq. Q can be quite 
conveniently carried out by inverse Laplace transform methods. The corresponding formulae 
can be easily read off from tables and have been collected for convenience in Ref. [1] 
The NNLL expression for the n-th moment can then be cast into the form 



* NNLL = j^pfe { "> 2 "> + 2 ^ h ' "> c *> "> «*•<*• "> } • («> 



5 



where q Hj i arises from the integration over the non-relativistic current correlator involving 
the dominant effective S-wave current and g n ^ from the current correlator involving one 
insertion of the 0(v 2 ) ~ (9(l/n)-suppressed S-wave current. Both correlators are referred 
to as Ai and A2, respectively, in Refs. [211 E2]. The contributions to g n> i coming from 
the Coulomb potential are UV-finite and therefore identical to the results given in Ref. [H] 
up to the previously mentioned trivial modification due to the Wilson coefficient of the 
Coulomb potential which differs from the strong coupling at NNLL order [32TI31] . The other 
contributions to g n) i as well as the result for Q n ^-, which contribute exclusively at NNLL 
order, differ due to different conventions for the potentials used in Refs. [27J E2] and also 
because they involve UV-divergences renormalized in the MS scheme (compared to a cutoff 
scheme used in Ref. [H]). All expressions for the non-Coulomb terms in g n i and for g n ^ 2 
have not been published before and are given explicitly in App. |Aj Throughout this paper 
all couplings and Wilson coefficients are understood to be renormalized in the MS scheme. 

The terms c% and C2 are the Wilson coefficients of the dominant and the 0(v 2 ) subleading 
currents, respectively. The NNLL order anomalous dimension of c\ contains the previously 
mentioned mixing and non- mixing corrections according to Refs. [2H - T26] and Ref. [27], re- 
spectively. These two types of NNLL order corrections to c\ have not yet been analyzed 
together in any phenomenological analysis of bottom pair production close to the threshold. 
A brief numerical analysis of C\ focusing on the NNLL order scale variation and the uncer- 
tainty due to the yet unknown soft mixing contributions in the anomalous dimension of C\ is 
carried out at the end of this section. Details on the analytic results for c\ and c 2 are given 
at the end of this section and in App. [B} 

We note that in our analysis we treat the charm quark as massless. Up to now the effects 
of the non-zero charm mass have only been treated in the fixed-order formalism exhibiting 
a shift in the bottom mass between —20 and —30 MeV [HI [35]. The corrections from the 
non-zero charm mass m c are conceptually interesting since m c is numerically close to the soft 
scale ~ a s rrib and generically above the ultrasoft scale ~ a 2 mb. The effects are therefore not 
suppressed by additional factors of m 2 /ml and are even somewhat enhanced since the charm 
mass interferes as a physical infrared scale in the renormalon cancellation. At this time a 
systematic renormalization group improved treatment of massive virtual quark thresholds 
in heavy quarkonium production is still lacking. However, taking the known effects in the 
fixed-order approach as a guideline, we conclude that the associated uncertainty is smaller 



than the perturbative uncertainties of our present analysis in Sec. IV 



The correlators and Wilson coefficients all depend on the dimensionless vNRQCD velocity 
renormalization evolution parameter v that is used to implement the correlated evolution of 
soft (/is) and ultrasoft scales (nu) according to the non-relativistic scaling \i v oc n 2 s /m h . For 
the large-n moments the typical choice for v is of order of the velocity of the bottom quarks 
which generically scales as \j^fn. The residual dependence of predictions on changes of v 
are used to estimate the perturbative uncertainty of the low-energy contributions contained 
in the theoretical prediction. On the other hand, changes coming from the variation of 
the matching scale /i^ can be used as an estimate of the perturbative uncertainty of the 
high-energy contributions contained in the moments. Both types of uncertainties can be 
considered as independent, because hard and low-energy contributions are separated in the 
non-relativistic effective theory. In order to avoid that the ultrasoft scale can exceed the 
soft scale we use the assignment \i$ — v^h, — z/2 A i h and impose the restriction v < 1. 
Moreover we parametrize the matching scale as fih — hm^, where the typical choice for h is 
of order one. The renormalization scaling parameter h has also been indicated in Eq. ([5]). 
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In our parametrization with the dimensionless renormalization scaling parameters v and 
h all explicit dependence on the bottom quark mass m& cancels from the Wilson coefficients 
as well as from the matrix elements corrections contained in g n ^ and q n ^- The only explicit 
dependence on the bottom quark mass arises from the overall dimensional factor shown in 
Eq. (|5]). In addition, there remains an implicit dependence on the bottom mass through 
the dependence of the Wilson coefficients or the matrix element corrections on the strong 
coupling a s , which has to be evaluated either at the matching, the soft or the ultrasoft 
renormalization scales. Because the matching scaling parameter h always occurs together 
with an additional m& factor and because the solutions of the anomalous dimensions can be 
expressed in terms of the strong coupling a s evaluated at the hard, soft or ultrasoft scales, 
it is straightforward to implement the /i-dependence of the moments by simply inserting a 
factor of h to the argument of each a s that arises. This method covers all /i-dependence 
except for two terms, one arising from the term ln(/i|/mf,/i[/) = In h in the NNLL non-mixing 
anomalous dimension of c\ [27] and one in the NNLL order matching condition ci( / u/ l , 1). 
Both terms are given in App. [Bias well. 

The moments shown in Eq. (pi have been written in terms of the bottom quark pole mass 
M^ ole , which is known to be sensitive to renormalon long-distance effects that seriously 
deteriorate the convergence of the perturbative expansion. It is therefore mandatory to 
employ a better defined short-distance bottom quark mass scheme. Short-distance mass 
schemes have an intrinsic dependence on an infrared cutoff scale [36], and for problems 
involving non-relativistic heavy quarkonium dynamics one has to employ so-called threshold 
masses, where this scale is of order of the inverse Bohr radius, 0{m\ ) a s ) [32]. In this work 
we will use the IS mass defined as half the perturbative contribution to the 3 Si bottonium 
ground state, T(IS') [37J|38]. Other possible threshold mass definitions have been discussed 
e.g. in Ref. [12], and in Ref. [22J the so-called renormalon-subtracted mass scheme was 
employed. The result can finally be converted to the MS mass, which is frequently used 
for predicting physical quantities involving scales larger than the bottom mass. Although 
this conversion avoids the renormalon problem of the pole mass, we note, however, that 
the conversion introduces a rather large dependence on the uncertainty in the value of the 
strong coupling a s due to a term cx a s rrib appearing in the conversion series. This motivates 
quoting the result for the IS mass as the main result of our analysis^] We note, however, 
that in our case this a s -dependence just compensates for the a s -dependence of the IS mass 
we determine from our sum rule analysis. 

To convert Eq. ^ to the IS bottom mass scheme its relation to the pole mass is required, 
which at NNLL order reads 

M b pole = Ml s {l + A LL + A NLL + [(A LL ) 2 + A™ LL + A^ NLL ]} . (6) 



The various A terms are given in Ref. [31], where it is necessary to use Eq. (A3) for the 

( s \ — i s \ 

coupling parameter a and the result for V^g ff from Ref. [32] for the coefficient . The 

logarithms in the A terms are to be understood as L = ln[/zs/(amf,)], which reduces to 

ln[z//a] for h = 1 as given in Ref. [31J . On the RHS of Eq. ^ the terms A LL and A NLL are 

the LL and NLL order contributions in the non-relativistic expansion, and the three terms 

in the brackets are the NNLL order corrections. As discussed in detail in Ref. [T7J, it is 

crucial to reexpand the series in a way that is consistent with the exponential non-relativistic 



1 We refer the reader to Ref. |39) for a review on the in our view incoherent situation concerning the value 
and uncertainty in a s (M z ). 
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expansion given in Eq. Q where the LL order contribution to the binding energy is treated 
as a leading order term. This is achieved by using the expression 



(M fc pole ) 2n (M 6 15 ) 2n 



' exp(-2nA LL ) 



x 



{l - 2nA NLL + n[(A LL ) 2 - 2A NNLL + 2n(A NLL ) 2 ] +...}, (7) 



where the first, second and third term in the second line correspond to LL, NLL and NNLL 
order, respectively. In order to consistently remove the infrared renormalon problem asso- 
ciated with the pole mass when switching to the IS scheme, it is important to consistently 
expand out the series in the curly brackets with the corresponding LL, NLL and NNLL order 
terms contained in p Uj i, while the exponential term containing A LL remains unexpanded. 
For the p„ )2 only the LL exponential term in Eq. (|7| has to be employed. Note that it is 
also crucial to use the same values for the renormalization scale parameters u and h in the 
Wilson coefficients and couplings of Eqs. ^ and ([7]). Unless stated explicitly we always 
refer to the IS mass definition as the bottom quark mass: = M fe is . 

The contributions to the renormalization group evolution of ci, the Wilson coefficient of 
the leading bottom pair production current, can be parametrized as 



ci(h, v) 

= 1 + e NLL + 



v) + (e NLL (^ ") + £„T L (^ "))+...] (8) 



IfNLlA 2 , tNNLL i tNNLL 
rjS J ' Sm ' ?] 



Mini 



+ . . . 



where £ NLL refers to the NLL order contribution and the £m NLL anc ^ (™ L to the NNLL 
order mixing and non-mixing corrections, respectively. The matching condition C\(h, 1) and 
the expressions for the £'s can be found in App. [B] and the references cited there. In the 
second line of Eq. ^ we have displayed the consistently expanded form where the terms in 
the squared brackets represent the NNLL corrections and the explicit dependence of the £'s 
on the scale parameters has been suppressed. We note that since the LL order anomalous 
dimension for the current is zero, the expanded form does not contain any exponentiated 
LL order contribution. 

Based on the recently completed calculation of the ultrasoft NLL running of the Wilson 
coefficients associated to the 0{v 2 )- and 0(a s v )-suppressed potentials [23H2S] , the ultra- 
soft mixing contributions to £ m NLL , referred to as Cm^oft nas Deen determined in Ref. |25j. 
Concerning the corresponding soft mixing contributions currently only those coming from 
the NLL order anomalous dimension of the spin-dependent (9(w 2 )-suppressed potential are 
fully known and were found to be tiny [28] . On the other hand, because the NLL matching 
conditions of all the suppressed potentials are known, it is possible to compute the term 
oc a^lnv of the soft mixing contributions, which we call £ m ^ofti- The result was already 



presented in Ref. [27] and it also given in Eq. (B7). 

The NNLL order non-mixing contributions in £nm LL ' i- e - the ultrasoft and soft contri- 
butions, are both completely known already since some time from Ref. [27] ■ In the same 
publication it was observed that the ultrasoft non-mixing contributions were more than an 
order of magnitude larger than the soft ones, and that the smallness of the latter was not 
arising from any accidental cancellation between different color factors but was a genuine 
property of all soft non-mixing contributions. In the following we analyze all known NNLL 
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contributions together numerically and give arguments suggesting that the ultrasoft con- 
tributions also dominate the mixing corrections and that the unknown soft corrections are 
negligible. 




FIG. 1. Panel a): RG evolution of the 3 Si current coefficient ci(y) normalized to c\{l) for rrn, = 
4.7 GeV and h = 1. The dashed line represents the full NLL result in expanded form: 1 + £ NLL . 
The solid line accounts for all known contributions including the new ultrasoft NNLL mixing 



correction in Eq. (B2) as well as the leading soft mixing logarithm at NNLL ^NIX 

in Eq. dB7|: 



1 + £ NLL + (£ NLL ) 2 /2 + ^nll + £NNLL ft + fNNLL^ The gray area aroun d the black line is generated 
by varying the known soft NNLL contributions to that curve by a factor between and 2. Panel b): 
Separate curves for the different (soft /ultrasoft, mixing/nonmixing) NNLL corrections (£ NNLL ) to 
the running of c\ as indicated in the plot. The NNLL contribution coming from the exponentiation 
of the NLL term, (£ NLL ) 2 /2, is tiny. Since it virtually coincides with the horizontal zero-line in 
this plot, it has not been drawn. The curve for Cnmsofti represents the linear logarithmic term 
ex aflnz^ contained in the complete soft non- mixing correction ^nmsoff ^ or both plots we have 
used oi n/=4) (4.7 GeV) = 0.217. 

In the left panel of Fig. [Tjthe renormalization scale evolution of c\(h, u)/ci(l, 1) is dis- 
played for 0.26 < v < 0.8 at NLL order (dashed line) and including all known NNLL order 
contributions, where the soft mixing corrections are approximated by Cm^ofti ( s °hd line). 
The gray band drawn around the NNLL result arises from the variation of the result coming 
from multiplying the soft terms at NNLL order by a factor between zero and two. Com- 
paring the width of the gray band with the overall shift between the NLL and NNLL order 
results and also with the overall remaining v dependence of c\(h, v)/c\(l, 1) at NNLL order 
we make two important observations. 

First, the ultrasoft contributions dominate the NNLL order result and are about an order 
of magnitude larger than the known soft contributions. From the right panel of Fig. [T] we 
see that this is true for the mixing and non-mixing contributions separately as well as for 
the sum. From the close vicinity of the curves for the full expression for soft an d its 



linear logarithmic contribution C^softi ^ a h^ ni> , given in Eq. (B9), we can also conclude 
that the curve for £m^ofti should represent the correct order of magnitude for the full set of 
soft mixing corrections. It is therefore reasonable to take the gray band as an uncertainty 
associated to the currently unknown soft mixing corrections. Comparing the gray band to 
the size and scale dependence of the bigger ultrasoft corrections we can conclude that this 
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uncertainty is smaller than the remaining NNLL renormalization scale variation we observe 
in the relevant range 0.3 < v < 0.6. This suggests that the unknown soft corrections can be 
neglected for a NNLL order prediction. 

As the second observation, we see that the NNLL order corrections to the anomalous 
dimension of C\ are about a factor two larger than the NLL order ones. This somewhat 
unsettling fact might cast slight doubts concerning the quality of the perturbative expansion 
due to large ultrasoft corrections, but might be as well associated to an anomalously small 
NLL order correction, which amounts to less than 10%. Nevertheless, given the situation it 
is certainly appropriate to invest a closer look on the behavior of the perturbation series for 
the moments. 

Finally we comment on the presence of nonperturbative power corrections in the OPE 
of the current correlators Ai- The effect of the leading power correction on the theoretical 
large- n moments in Eq. ^ is associated with the (9(Aq CD ) gluon condensate arising in the 
OPE for A\. The corresponding nonperturbative power correction to the n-th moment can 
be written approximately as [7j 



n 3 7r 



8p th,n P = pth,LL {ag G 2 } « exp( _ Q 4 ; (9) 

where we take a s in the exponent at the soft scale. To our knowledge the subleading correc- 
tions to the Wilson coefficient of the gluon condensate are currently unknown. Numerically 
the relative correction from Eq. d9|) is less than a percent for n < 20, i.e. even smaller then 



the experimental uncertainty and therefore negligible for our purpose, cf. Sec. IV In order 
to make sure that also higher order power corrections [ID] are sufficiently suppressed, we, 
however, only consider n < 10 for our analysis. 



III. THE PERTURBATIVE EXPANSION 

The perturbative series for the large-n moments P* h in Eq. ^ follows the non-relativistic 
expansion scheme of Eq. ^ and, as explained in the previous section, requires that the 
perturbative series for p n l coming from the leading order current correlator is consistently 
expanded with the additional corrections of Eq. (|7j) that arise from employing the IS thresh- 
old mass scheme. However, concerning the factor c\ in Eq. ^ one has the option either to 
keep it as an overall factor for the series associated to Q n> i ( "factorized form" ) or to expand 
it out as well together with the series of Q n) \ ("expanded form"). The difference between the 
two ways to expand P* h represents terms from beyond NNLL order and should, in principle, 
be within the theoretical uncertainty of our NNLL result. However, given the observations 
concerning the perturbative behavior of c\ discussed above we believe that it is mandatory 
to have a closer look on the expanded and factorized expansions of the moments. 

An important property of the large-n moments is that the non-relativistic dynamics 
encoded in them is governed by bottom quark velocities v ~ 0(l/^/n). Since the vNRQCD 
renormalization scaling parameter v is typically of order v , we would expect that a good 
convergence of the series and the best physical description is achieved for v ~ (l/y/n + 
const.), where the additional constant term arises from a saturation effect that protects the 
renormalization scale from vanishing for large Next we examine for the expanded and 
the factorized form of the large-n moments, how well they satisfy this behavior. 

2 As argued in Ref. [7] the proper choice for the soft {us) and ultrasoft renormalization scale (fijj) m the 
limit n — > oo approaches the inverse Bohr radius (~ a s mf,) and the binding energy (~ a 2 s rrib) of the lowest 
lying resonance, respectively. 
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In Fig. [2]we plot the values of the renormalization parameter v that are required to satisfy 
v) = P% xp as a function of n using h = 1 and a fixed value of m^, where the P£ xp are 
the experimental moments. The value of the bottom quark mass has been fixed by fitting 
the theoretical 10-th moment, P*q> t° the experimental one using v = /(l/vlO + 0.2). Here 
the constant / is a fudge parameter which should be of order one and which we have varied 
between 0.8 and 1.25 in steps of 0.05. The results for v we obtained for n ^ 10 are shown as 
blue dots, and values for equal / are connected by the blue lines. The left panel shows the 
outcome for the expanded moments and the right panel for the factorized moments. We see 
that the v values for the fully expanded moments decrease with n and are indeed consistent 
with the expected behavior v ~ (l/y/n+ const.). On the other hand, the v values from the 
factorized moments increase with n and are inconsistent with the expected behavior. As a 



comparison we have displayed the functions v 
dashed red lines. 



/(l/Vn + 0.2) with / = 0.8, 1.0, 1.25 as the 




0.3 L 



10 
n 



12 



14 




FIG. 2. Values for the renormalization parameter v required by P^ h (mi,; u) = P% xp for fixed mass 
nib and different n. The blue lines connect the v values (blue dots) for the same reference mass. 
Lines for 10 different reference masses determined from P$(m b ; v = /(l/\Zl0 + 0.2)) = P^ p with 
/ between 0.8 and 1.25 are shown. For the plot in panel a) we used the fully expanded NNLL 
expression and for the plot in panel b) the factorized NNLL expression for P^ h , as described in the 
text. The reference masses in a) range from 4.77 GeV for / = 0.8 to 4.74 GeV for / = 1.25. The 
reference masses in b) range from 4.69 GeV for / = 0.8 to 4.64 GeV for / = 1.25. The dashed red 
curves are defined by v = f(\/y/n + 0.2) with / = 0.8, 1, 1.25 respectively. 

The quite different outcome for the expanded and factorized moments is caused by the 
perturbative behaviour of c\ discussed in the previous section. This shows that the (ex- 
pansion scheme dependent) fraction of terms from beyond the NNLL level included in the 
two ways of the NNLL expansion for the moments, are large enough to qualitatively affect 
renormalization scaling features of the perturbative result. Because only the fully expanded 
moments show the expected behavior for their renormalization scaling parameter u, we 
adopt them for our numerical analysis to determine the bottom quark mass. However, we 
require as an important cross check that the bottom mass values determined from the fac- 
torized moments have to be consistent with those from the expanded moments within the 
perturbative uncertainty. 
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IV. NUMERICAL ANALYSIS 



The experimental moments we use for our analysis are based on the measurements of the 
masses and electromagnetic decay widths of the first four T resonances (T(15)-T(45)) [13] 
in the narrow width approximation, BABAR data for the threshold region [4T] in the en- 
ergy range between ^/s = 10.62 and ^/s = 11.21 and the perturbative QCD result for the 
continuum at energies ^fs > 11.21 |42j . where no experimental data is available]^] As we 
are aiming at an analysis for large-n moments, the T resonances and the BABAR region 
constitute the major part of the experimental moments ((87,93,96,98)% for for the T reso- 
nances, and (5.7,3.7,2.3,1.4)% for the BABAR region for n = (6,8, 10, 12)) and only little 
effort has to be invested for the continuum region above 11.21 GeV where no experimental 
data exists. We follow the approach of Ref. [13] for the BABAR region and also adopt an 
additional 10% model uncertainty for the contiuum region. We stress, however, that the 
latter uncertainty only constitutes a very small and numerically irrelevant part of the whole 
error budget of the large-n moments which are dominated by the uncertainties from the T 
resonances and the BABAR region. We note that applying our approach for low-n moments, 
we obtain results that are perfectly compatible with the moments given in Ref. [B] albeit 
with somewhat larger errors if applied for low n values, as in their approach only pertur- 
bative QCD uncertainties were accounted for in the continuum region. The mean values 
together with the respective statistical and systematic errors of the P% xp for 4 < n < 14 are 
given in App. [CJ 

In order to determine the IS bottom mass from single moment fits we simply solve the 
equation P^ h (mb) = P% xp for m b . In order to estimate the perturbative error we carry out 
these fits many times varying matching and renormalization scales in reasonable ranges. To 
parametrize the scale variations for the different moments P^ h such that the n-dependence 
of the renormalization scales and their correlation is implemented coherently we use the 
assignments 

/j, h = hm b ,fi s = hm b fv* and \l v = hm b (fv*) 2 (10) 

for the matching, soft and ultrasoft scales, respectively, where v* = (l/y/n + 0.2) is our de- 
fault choice for the renormalization scaling parameter v. So the matching, soft and ultrasoft 
scales adopt their default values for h = f = 1, and we parametrize variations of these scales 
by scanning over a certain region in the two-dimensional h-f '-plane. 

The contour plot in Fig. [3] shows the results for the IS bottom mass obtained for n = 10 
at NNLL order using the expanded expression for the theoretical moment in the h-f -plane. 
The default values are h = f = 1 (red dot). For estimating the perturbative error we vary h 
and / around their default values applying the restriction that the ultrasoft renormalization 
scale remains within the range 0.5 m^z/* 2 < flu < 2m b v* 2 . Furthermore we fix the range for 
the variation of h to 0.75 < h < 1/0.75. The two conditions define the region indicated by 
the red dashed lines, and we note that for both h and / being small this region restricts the 
ultrasoft scale from dropping below 0.13 m b for n = 10. 

Figure [3] displays the overall behavior for the bottom mass extracted from the NNLL order 
moments, indicating steeply rising mass values for decreasing values of h and / corresponding 
to the matching, soft and ultrasoft scales being small. However, in the restricted region we 



3 We thank V. Mateu and B. Dehnadi for providing us the data compilation and a numerical code for the 
experimental moments. 
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FIG. 3. Contour plot of the IS bottom mass determined from P^q (m&) = P\q P as a function of the 
parameters h and / as defined in the text. The different contours are labeled by the respective mass 
value in GeV. The region in the h-f plane bounded by the red dashed line represents the parameter 
space we scan to determine the variation of the mass, which contributes to our perturbative error 
estimate. The region is defined by 0.75 < h < 1/0.75 and demanding that 0.5 fj^j < /J,u < 
where \x{j = nribv* 2 ■ The red point inside this area indicates our default values / = h = 1 for the 
mass determination. 



use for estimating the perturbative error, the contours are relatively flat with the default 
rrib representing a central value within the overall range of values obtained in the restricted 
region. 

For comparison we have also displayed in Fig. [4] the corresponding bottom mass contour 
plots obtained from the LL order moment (left panel) and from the NLL order moment 
(right panel). It is a conspicuous fact that in the restricted h-f region the form of the 
contours and the range of the bottom masses covered for the LL analysis is quite similar to 
the NNLL order analysis. On the other hand, for the NLL order analysis the outcome is 
quite different, since the mass values strongly decrease in the small h-f region and since the 
default bottom mass is very close to the upper end of the obtained range of bottom mass 
values. Moreover there is hardly any overlap to the mass range obtained from the NNLL 
and LL analyses. The overall behavior we see from the outcome of this analysis at LL, NLL 
and NNLL order confirms the slightly marginal character of the perturbative series for the 



large-n moments we already discussed in Sec. Ill 



From our observations it is clear that one has to conclude that the LL and NLL or- 
der renormalization scale variations do not provide reliable estimations of the perturbative 
uncertainty at their respective orders and that the uncertainties are in fact substantially 
larger. On the other hand, the consistency of the analysis at LL and NNLL order and 
between the NNLL order bottom mass results for the expanded and factorized moments, 
as shown below, are indications that the NNLL analysis is reliable and that the variations 
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h h 

FIG. 4. Contour plots of the IS bottom mass contour plots determined from the 10th moment 
at LL (left panel) and NLL (right panel) order. The gray region in the lower left corner of the 
NLL plot indicates that the equation P 10 ' (m&) = P^ p does not have a solution related to the 
ultrasoft scale being too small (flu < 0.6 GeV). 

of / and h in the region described above give a proper estimate of the perturbative NNLL 
order uncertainty. We therefore adopt the NNLL order large-n moments and the procedure 
as described above to determine the bottom mass and the perturbative uncertainty for our 
final numerical analysis following below. 

In Fig. [5] we present the results for the IS bottom mass from the moments for n = 
6,8,10,12 in the expanded form as a function of a s (Mz) at NNLL order (red) and for 
comparison also at LL (blue) and NLL order (orange). Given the discussion above we will, 
however, only discuss the NNLL order results. The respective dashed lines represent the 
mass for the default scale choices with f — h — 1 and the bands come from varying / 
and h in the region as discussed in detail above. The masses obtained from the default 
scale choices are remarkably consistent for the different n values exhibiting deviations of 
less than 20 MeV for a s (Mz) in the region around 0.118, where the current world average 
a s (Mz) = 0.1183 ± 0.0010 |39J, is located. The width of the band from the scale variations, 
on the other hand, slightly decreases for increasing n indicating that the dependence of the 
moment normalization on the exponent of the mass slightly overcompensates the increase 
of higher order a s corrections due to the smaller renormalization scales. Overall, the IS 
mass slightly increases with a s , but the dependence is rather weak and linear to a very 
good approximation so that the error in a s (Mz) only has a rather small impact on the 
final uncertainty in the bottom mass. We note that the uncertainties in the bottom mass 
coming from the error of the experimental moments amount to (6.5,4.3,3.3,2.7) MeV for 
n = (6,8,10,12) and are much smaller than the perturbative error and the uncertainty 
coming from a s . For quoting the final result the experimental uncertainty does not play any 
essential role. 
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FIG. 5. IS bottom mass with errors determined from P*' l (mb) = P% ep for n = 6,8,10,12 as a 
function of a s {Mz) using the strictly expanded expression for P^ h (mb). As indicated in the panel 
for n = 6 we plotted for each moment three (partly overlapping) bands with LL (blue), NLL 
(orange) and NNLL (red) precision, respectively. The error bands are generated by varying the 
parameters / and h of the theoretical moments P^ 1 within the parameter space defined in Fig. [3] 
and adding the corresponding experimental error in quadrature. 



In Fig. [6] we show the corresponding results for the IS bottom mass using the factorized 
expression for the theoretical moments. Here the LL results are identical to the previous 
analysis based on the expanded moments. While the width of the NNLL order bands from 
scale variations are almost identical to the expanded moment analysis, we see that the NNLL 
order results for the mass values are somewhat smaller, but consistent within errors with 
the results for the expanded moments. The n-dependence of the mass values for a s values 
around the world average is about three times larger than for the fully expanded moments. 
We note again, that the factorized moments have an inconsistent renormalization scaling 
behavior, and we therefore adopt the moments in the expanded form for the main numerical 
analysis. The consistency of the results from expanded and the factorized moments, on 
the other hand, indicates that the two results are in agreement within their perturbative 
uncertainties. 

For quoting our final number for the IS bottom quark mass we take the default value 
obtained from the 10-th NNLL order moment in expanded form as the central value and 
adopt half the size of the band from scale variation as the perturbative uncertainty. Using 
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FIG. 6. Band plots analogous to Fig. [5] but with the factorized expression (as defined in the text) 
instead of the expanded expression used for the theoretical moments. 



the current world average for the strong coupling a s (Mz) = 0.1183 ± 0.0010 [39] we obtain 

Ml s = 4.755 ± 0.057 pcrt ± 0.009 Qa ± 0.003 cxp GeV. (11) 

It is worth to emphasize that the a s error increases by 9.3 MeV for each unit of 0.001 the 
uncertainty in a s {Mz) is increased. It is also interesting to quote the result as a function 
of a s {Mz)- Since the dependence of the bottom mass on the value of the strong coupling 
is linear to a very good approximation, we can quote the a s -dependence of the bottom IS 
mass by the following numerical fit formula: 



M, 



is 



4.755 + 



a s {M z ) - 0.1183 
0.107 



± 



0.057 + 



a s (M z ) - 0.1183 
037 



pert 



GeV, 



(12) 



where we have dropped the experimental error. The deviation of the linear approximation 



in Eq. (jT^J) to the exact result is less than 0.6 MeV for 0.113 < a s (Mz) < 0.120. 

It is straightforward to convert our result for the IS mass to the MS mass by combining 
the pole- IS mass relation in Eq. (|6 



and the pole-MS mass relation. For this conversion 
one has to employ the e-expansion [TTl 1371 138] to consistently cancel the renormalon con- 
tributions. The systematic cancellation of the renormalon contributions also requires that, 
within the conversion, the corrections in the pole-lS mass relation are evaluated at the same 
renormalization scale as those in the pole-MS mass series. This can be achieved with the 
best convergence using the hard scale M 6 1S and requires to use the fixed order expansion 
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for the pole-lS mass relation which is obtained from Eq. ^ for v = 1 and h = 1. Since 
the conversion formula has a significant a s -dependence due to a term oc a s m b arising in the 
series relating the MS and pole masses one needs to convert using the full a s -dependence of 



the IS mass result as displayed in Eq. (12). 



Using again the world average for the strong coupling as an input we then obtain 

m b (m b ) = 4.235 ± 0.055 pcrt ± 0.003 cxp GeV, (13) 

where we have added an additional conversion error of 15 MeV to the perturbative uncer- 
tainty coming from half the size of the three-loop correctiorj^jin the perturbative series. The 
resulting value for the perturbative uncertainty is, however, still smaller than for the IS mass 
due to the conversion series. The uncertainty coming from the value of a s reduces to less 
than 0.5 MeV, as the sizable a s dependence in the 1S-MS conversion series mentioned above 
is anticorrelated to the monotonically increasing a s dependence of the IS mass. Remarkably 
these two effects cancel almost entirely. We have therefore dropped the a s induced error 



from Eq. (13). In analogy to Eq. (12) we can also give the result for the MS mass showing 
the full a s dependence linear fit function: 

a s (M z ) - 0.1183 



m b (m b ) = 4.235 ± 



0.055 + 



GeV, (14) 

pert 



0.41 

where we did not include a linear a s fit term for the central value as it contributes less than 
one MeV for 0.113 < a s (Mz) < 0.120. Again we have dropped the experimental error. 

Our result for the bottom MS mass m b (m b ) should be compared with the corresponding 
result obtained in the analysis of Ref. [22]. Their analysis was also based on RGI large-n 
moments having NNLL order input, but was using only the NLL order result for the anoma- 
lous dimension of the leading order current Wilson coefficient c\. Thus their theoretical 
moments were missing the numerically important NNLL contribution to the evolution of 
c\ shown in Fig. [T] Since c\ enters as a square, their theoretical moments were by about 
30% lower than ours. For a single fit for the 10-th moment one therefore expects, due to 
the l/m^ n overall dependence on the mass as shown in Eq. ([H]), that their mass result is 
by about 1.3% lower than ours. This amounts to about 50 MeV and is consistent with the 
result m b (fn b ) = 4.19 ± 0.06 MeV obtained in their analysis. The perturbative uncertainty 
estimated in their work was quoted between 45 and 65 MeV and is numerically similar to 
ours. We note, however, that they determined the uncertainty from the difference in the 
bottom mass obtained from the NLL order and their incomplete NNLL order moments. 

We also note that our result for the MS mass is not quite compatible with the results 
given in Ref. [45] (jn b (jn b ) = 4.163±16 MeV) or even in Ref. [TO] (jn b (jn b ) = 4.171±9 MeV) 
which were obtained from 0(a^) fixed order analyses related to low-n T sum rules and 
taking their uncertainties exactly as quoted. Here we note again that we have treated the 
charm as a massless quark. In previous fixed-order analyses [THl E5] it was found that the 
non-zero charm mass causes shifts in the bottom quark mass obtained from large-n sum 
rules between —20 and —30 MeV. This sizable corrections is related to the fact that the 
charm mass is similar to the inverse Bohr radius a s m b and therefore not l/m^-suppressed. 
At this time the renormalization group improved non-relativistic effective theory approach 
has not yet been extended systematically to account for massive virtual quark thresholds. 
The required conceptual developments are theoretically interesting and might contribute 
reconciling the discrepancy. For comparing our result to bottom mass determinations from 
other physical quantities we refer to Ref. 



This corresponds to the 0(e 3 ) terms in the epsilon expansion, see e.g. Eq. (99) of Ref. 
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FIG. 7. Comparison of the masses obtained from the fixed order (a) and RGI calculation (b) of the 
10-th moment, P*q (m&) = -Pfo^- ^ n P ane l s a anci b we show the mass values with LO, NLO, NNLO 
and LL, NLL, NNLL accuracy, respectively. The corresponding error bands were generated in the 
same way as in Fig. [5] Concerning panel a, we note that for some low raj values in the NLO band 
and the associated values for h and / the ultrasoft coupling as(fJ-u) reaches 0.65 causing numerical 
instabilities. 

V. FIXED-ORDER MOMENTS AND MULTIPLE MOMENT FITS 

Before RGI predictions based on the extended versions of NRQCD were possible a num- 
ber of fixed-order analyses were carried out at the NNLO level [HHUj. As mentioned in 
Sec. [T| the fixed-order moments showed a rather bad perturbative behavior, and the different 
analyses addressed this issue in different ways in order to achieve bottom mass determina- 
tions with small uncertainties. We can turn our RGI moments into fixed-order moments 
by setting the matching, soft and ultrasoft scales equal to the soft scale except for the 
current Wilson coefficients, where we keep the matching scale fih- This leads to a logarithm 
ln(/i5/7//j) = ln(/i/*) appearing in the current Wilson coefficient c\ whose coefficient does not 
depend on the regularization scheme nor the conventions used for the potentials. Although 
the form of the resulting moments do not agree exactly with either one of the analyses in 
Refs. [T4T - fT8"] due to differences in the regularization scheme, the expansion scheme or in the 
conventions for the potentials, the resulting fixed-order moments provide a good insight into 
the improvement obtained by the renormalization group logarithms summed into the RGI 
moments. 

The IS bottom mass results obtained by carrying out a single moment fit for our fixed- 
order moment for n = 10 in the expanded form and using the same variations of h and / we 
employed in the RGI analysis^ which now only cause correlated variations in the matching 
and the soft renormalization scale, is shown in Fig. [7j The result should be compared with 
the outcome using the RGI NNLL order moment P{q shown in the lower left panel of Fig. |5j 
which we have for convenience displayed once more in the right panel of Fig. [7} Because 
at LL order the soft scale is the only renormalization scale that arises in the theoretical 
moment, the LL RGI and LO fixed-order results agree exactly. Substantial differences are, 

5 The region in the h — f plane we use in our analysis leads to the variations 0.75m& < fi^ < 1.33m;, and 
0.32m;, < /is < 0.84mb for the matching and soft scales, respectively. For the matching scale the variation 
is somewhat smaller than what was employed in the fixed-order analysis of Refs. |14H18j but roughly 
agrees with those analyses for the soft scale. 
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however, visible for the NLL vs. NLO results where the fixed order result shows even lower 
bottom mass values and a substantially larger scale variation. Comparing the NNLL and 
NNLO results we see that the fixed-order results cover larger bottom masses and exhibit 
scale variations that are twice as big as for the RGI results. We clearly see the improvement 
related to summing the logarithmic contributions within the RGI approach. Determining 
the final result for the IS bottom mass from the NNLO moment in the same way as we 
did it in our NNLL RGI analysis we would obtain for a s (Mz) = 0.1183 ± 0.001 the result 
M b is = 4.791 ± 0.097 pcrt ± 0.012^ s ± 0.003 exp GeV (and thus fn b (m b ) = 4.269 ± 0.090 P e rt ± 
0.002 Qs ±0.003 exp GeV for the MS mass), which is compatible with our results from Eqs. (TTj) 
and (13), but with a substantially larger perturbative uncertainty. Our fixed-order result is 
in agreement with Ref. p2] (m b (fn b ) = 4.26±0.09 pert ±0.01 Q , s ±0.02 exp GeV) which also used 
a NNLO fixed-order single moment fit for n = 10 and at the time of the analysis accounted 
for larger strong coupling and experimental errors. 

Another instructive analysis is related to carrying out simultaneous fits to several mo- 
ments using a x 2 analysis. It is straightforward to carry out such fits with our moments 
as well. In Fig. [8] we show the allowed bottom mass ranges obtained from x 2- fits based on 
the moments for n = 4,6,8, 10 where brown refers to the fixed-order NNLO analysis and, 
for comparison, red to the RGI NNLL order analysis. The upper panels show the outcome 
for the n-dependent assignment of the renormalization scales and the variation in the h — f 
plane as described for the previous analyses, i.e. the choice of scales for the moments entering 
the x 2 function is n-dependent. Here the left panel is obtained using the full experimental 
covariance matrix (based on treating the experimental input data for masses, widths and 
i?-ratio values as statistically independent), and the right panel is obtained treating the 
moments as statistically independent. The result for independent moments agrees almost 
exactly with the result of the n — 10 single moment fit, which is easily understood since 
Piq P has the smallest error and thus the highest weight for uncorrelated moments. 

Accounting for the correlation, on the other hand, leads to a different outcome, where 
the NNLL order bottom mass values are lower and the NNLO fixed order results lead to 
much smaller scale variations compared to the single moment or the uncorrelated fits. We 
have traced this behaviour back to the fact that the correlation coefficients are very close to 
unity and that the x 2 function fits the n-dependent shape of the moments rather than the 
actual values of the moments. Together with the correlated scale dependence of the theory 
moments this leads to the observed strong cancellation of the scale variation. This is also 
associated to rather large minimal x 2 /dof values with an average of around 50. From this we 
conclude that the correlations are too large and the resulting scale variations are unreliable. 
Accounting in addition for the fact that the theory uncertainties are much bigger than the 
experimental ones, we conclude that single moment fits seem to be the only practical option 
to extract the bottom mass. 

This conclusion is also reconfirmed by the analysis shown in the lower two panels. Here 
we have displayed the results for the bottom mass using the same x 2 fit method as for the 
respective upper two panels, but now using the renormalization scales adopted for n = 7 (the 
average of the n values used for the fit) for all moments. Treating the experimental moments 
as uncorrelated (right lower panel) the results for the bottom mass are very similar as for 
the moments with n-dependent scales (shown in the right upper panel) since again the 10-th 
moment has the highest weight in the fit. On the other hand, including the correlation we find 
that almost all scale variations cancel. For the NNLL moments the resulting scale variation 
is between 25 and 30 MeV, while for the NNLO moments it is below 5 MeV. This remarkable 
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FIG. 8. Multiple moment fits for the IS mass with n = {4, 6, 8, 10}. Each panel shows two bands 
associated with a fit using the NNLO fixed order moments (brown) and the NNLL RGI moments 
(red), respectively. The bands are generated by scanning the usual h — f region (see Fig. [3]). Fits 
scanning an n-dependent v range with v* = (l/y/n-\- 0.2) and an n-independent v range with 
v* = (l/y/7 + 0.2) while allowing for full/no correlation of the experimental moments are shown 
as indicated by the labels of the individual panels. 



cancellation can be traced back to the fact that the experimental correlation matrix puts 
the highest weight to a linear combination of moments that minimized the sensitivity to the 
uncertainties of the T electronic widths, which represent the biggest source of uncertainty 
for the large-n experimental moments. Using the same renormalization scale values for all 
moments reduces very effectively the sensitivity to the electronic width of the perturbative 
T bound states entering the theoretical moments. Since the theoretical electronic widths, 
and the theoretical i?-ratio in general, are the biggest source of scale variation in the theory 
moments, a common choice of renormalization scale in all moments, leads to the observed 
reduction of scale variation. 

Given our observation that the experimental correlation of the moments is too strong to 
allow for reliable multiple moment fits, we have to conclude that the tiny scale variation 
shown in the lower left panel is an artifact of the scale choice and does not represent an 
estimate of the perturbative uncertainty. Overall the unreliable multiple moment fits are a 
consequence of the fact that the theory errors are significantly larger than the experimental 
ones. We note that fixed-order multiple moment fits with n-independent renormalization 
scales based on several different sets of moments were used in Refs. [T4"t \T7\ |T8] . The 
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theoretical uncertainties quoted in these analyses originated mainly from using the different 
sets of moments and not from the scale variation. 



VI. CONCLUSION 

In this work we have determined the IS bottom quark masses from large-n T sum rules 
with renormalization group improvement. Since for large values of n the quark dynamics en- 
tering the moments is non-relativistic, the experimental moments are dominated by the well 
known data on the T resonances and the experimentally unmeasured continuum contribu- 
tions are not relevant. The theoretical computation of the moments involves a simultaneous 
expansion in a s and 1 / y/n and is carried out within the framework of non-relativistic effec- 
tive field theories. In general, within the non-relativistic framework where n is large, one has 
to treat terms ex a s / ' \Jn and oc a s ln(n) of order unity. For the renormalization group im- 
proved treatment both types of terms are summed to all orders following the generic scheme 
of Eq. (|3]). We have carried out this scheme using the vNRQCD ("velocity" NRQCD) for- 
malism originally devised in Ref. j2TJ and extended in a number of subsequent publications 
as discussed in the main body of this paper. We have also determined the MS bottom quark 
mass m^ijn^). We note that the conversion involves a sizable dependence on the value of 
the strong coupling a s . While we find that the bottom IS mass slightly increases with a s , 
the MS mass m^mi,) turns out to be essentially c^-independent. 

In this work we have for the first time carried out a renormalization group improved anal- 
ysis using all available NNLL order results. At the NNLL order level currently all corrections 
are known except for the so-called soft mixing corrections to the anomalous dimension of the 
leading order heavy quark pair production current. We have given arguments showing that 
the known full set of NNLL order ultrasoft corrections numerically dominate over the soft 
contributions and that the uncertainty coming from the set of unknown soft corrections are 
negligible compared to the large ultrasoft contributions and their remaining scale variation 
at the NNLL order level. Our results can therefore be considered as the outcome of a full 
NNLL order analysis. 

Prior to our work a number of large-n analyses were carried out, which were based on 
the so-called fixed-order approach, where only the terms oc a s / \fn were summed up system- 
atically [HHT7]. One important characteristic of these fixed-order analyses was that they 
exhibited rather bad convergence properties. We found that the summation of the loga- 
rithms contained in our moments improves the convergence considerably. Nevertheless the 
perturbative series for the renormalization group improved moments shows sizable correc- 
tions which made a careful examination of the perturbative uncertainty of our NNLL order 
prediction mandatory. 



Our final results have been presented in Eqs. (11-14) of Sec. IV 
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Appendix A: NNLL calculation of the moments 



The calculation of the contributions g n> i and p„ )2 to the n-th NNLL moment in Eq. ^ 
follows the lines of Ref . [H] . The results can be written as 

Qn,l = [Qn,l]c + [f?„,l]kin + (yi S \v) + 2V«(")) [g ntl ) S + V«(l/) [ Qn ^] r 

c 2 

— CaCfOL s [fe,l]cACF + -TT®S [f?ra,l]cF2 



and 



+ <4 Viai") [Qn,l]kl + ®sV$(v) [Qn,l]k2 



2 ■»)(*)/ 



Qn.2 



An 



2+ -v^<p-4V^Y1 

6 P =2 



, p 2(p-3) C P 



r(¥) 



(Al) 



(A2) 



We abbreviate the strong coupling at the matching, soft and ultrasoft scales as ah = a s {Hh), 
a s — ®s(Hhfv*) and an = a s ([ihf 2 v* 2 ), respectively, and define 



a \/n 



a = 



1 

47T 



v c ( 5 ff (^) 



(A3) 



The expressions for the potential Wilson coefficients, V^, can be found in Ref. |32j. Except 
for the contribution of the Coulomb potential [g n ,i ] c , which contributes at LL and NLL 
order, all the other terms contained in Eqs. ( A1| A2) constitute NNLL order contributions. 
The LL Coulomb contribution reads 



1 + 2^0 + 4^ 



p=2 



(A4) 



and differs from the fixed order LO result given in Ref. [14] only through the definition of 
the variable cf) which in the fixed-order treatment involves the replacement a — > CFa s (fj,s)- 
The difference constitutes a NNLL order correction. For the analytic expressions of the 
NLL and NNLL contributions to [f? n ,i]e coming from the one- and two-loop corrections to 
the Coulomb potential we refer the reader to the expressions (57-63) in Ref. [H], where the 
identification /i so f t = fis is understood. 

The contribution to the Coulomb potential arising from the QED photon exchange 
is formally a NLL order correction, when we count the electromagnetic coupling a Em as 
0(a 2 s ). The correction can be i mple mented in a straightforward way with the replacement 
V c ]f — > Vpg ff — AttQI a em in Eq. (A4). However, numerically one has to compare Cpa s with 
Q b a c 



, and for the strong coupling at the soft scale one obtains that {QbOt cm ) / {Cpois) 2 is 
(9(1%) which makes QED effects entirely negligible for our NNLL order analysis. For an 
implementation of QED effects a fixed-order treatment is therefore well justified. Apart 
from the QED effect in the Coulomb potential one then only has to consider a NNLL order 



correction in the matching condition for the current coefficient C\ given in Eq. (Bll ) through 
the replacement CVa s — > Cpa s + QbQ. e m in the NLL order correction. We find that this af- 
fects the bottom quark mass by less then 1 MeV, which confirms the observation made in 
Ref. |22]. 
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The remaining non-Coulombic contributions to Q n) \ differ from Ref. [14J because we use 
the differing conventions for the potentials from Refs. [27J |32j and dimensional regularization 
as the regularization scheme. The corresponding results are new and read 



[On, 



1 kin 



n 



p=2 



(p-l)(p-3) C P 



' p— 1 1 



+ \ W) + ( t ~ ln2 + H»s/m b ) ) [g n ,i}c L , 
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with the function f n (4>) defined as 
/n(0) = 2 {-l + y +ln(2v^) + 2v^0 + ln 



(All) 



oo >. 

2^y 

p=3 1 V 2 / 
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We note that for the NN LL o r der co rrection terms in the moments it is sufficient to use 
a ~ a LL = C F c^ L in Eqs. (|A^-(|ATl]). 



23 



Appendix B: Current Wilson Coefficients 



The NLL anomalous dimension of the 3 Si current coefficient C\(y) is given by [52" } |4H ] |4"T] 

NLL 



d 



V ( C S) M 



+ 2^c» 



(Bl) 



with S 2 = 2 for the heavy quark pair produced in the spin triplet state. The solution at 
NLL order, £ NLL , has been given in Eq. (57) of Ref. [32]. The form of NNLL ultrasoft mixing 
corrections in Eq. (B2) below implies the convention that the z and u = 1/(2 — 2;) parameters 
in the expression for £ NLL , as given in Ref. [32], are defined according to Eq. (|B6[ ) below. 

The ultrasoft contribution to the NNLL mixing correction £^ NLL in Eq. ([8j) was deter- 
mined in Ref. [23] from Eq. (Bl) through the subleading NLL ultrasoft RG evolution of the 
potentials V2,V r and V^s 



(see also Ref. ESI) and reads 



Sm,usoft 03 h 

Po 



7 ^ 

- H Vz 1 

4 6 



In 



57T 



In 2 



3-2z 



, h In 2 , 
2/ V 2 

41n(2 - z) 



+ z< 



2Li, 



1 



In z 



with 



A 
B 

z = 



— Ca Trif , Bi 

3 3 n H 



34 



C\ - AC F Tn 



20 



-C F (C A + C F )(C A + 2C F ) 
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(B2) 

(B3) 
(B4) 
(B5) 
(B6) 



The complete result for the soft NNLL mixing correction, £msoft * s currently unknown. 
However, its linear logarithmic contribution oc a\ In v was determined in Ref. [27] from the 
known subleading matching conditions for the potentials [1H1 H9] that appear in Eq. (Bl) 
and reads 



tNNLL 
?m,softl 



a 



48tt 



h Cl 



C A (l6 S 2 - 3) + AC F (5 - 2 S 2 ) - y Tj 
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(B7) 



The NNLL non-mixing contributions (™ LL were determined in Ref. [27J for the assignment 
fiu = fj^/m b , i.e. for h — 1. For the more general case \iu = (^/(hrrib) we use in our work 
(see the discussion in Sec. |TT|) , the NNLL order non- mixing anomalous dimension receives 
an additional additive term oc In h which reads 

4CV {C\ + 3C A C F + 2C 2 F ) 



3/3 



a, 



ln(/t) (2 In 



1 

2 - z 



-z + 1). 



(B8) 
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For a comparison carried out in our numerical analysis we also give the linear soft non-mixing 
logarithm contained in (™ L for h—1, which is determined by expanding the full result for 

tNNLL • n . 
"Snm.soft 111 a h- 



tNNLL 
?nm,softl 



a. 



288vr 



303C 2 A + [8C f (S 2 -12)-111Ca]/3o + 2C A C F (U1 -40S 2 )+360C| 



In v . 

(B9) 



Taking £m^ofti as an approximation for the complete soft NNLL mixing corrections, we can 
thus write the NNLL contribution to ln[ci(z/)/ci(l)] for arbitrary h as 



NNLL 



rfNJNLL/T _ i\ | XtNNLL/T \] , r f NJNLL , /-N IN LL ] 
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(BIO) 



The matching condition for the current Wilson coefficient ci(l) is given in Ref. [27] for 
h = 1. For arbitrary h we have to add a ln(/t) term which arises from the term ln^^/mf,) 
and find 
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1 a h + 
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cl 



ln2 31 2 

IT ~ 24 ~ ^ 
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where the constant k was defined in Eq. (21) of Ref. 



(Bll) 



The coefficient of the t> 2 suppressed 3 Si current c%(y) is needed at the LL level and 
reads 131 1 



, N 8Cf , ,ocux 1 
3p oi h 6 
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Appendix C: Experimental moments 



The experimental moments we use in this work together with their statistical and sys- 
tematic errors, respectively, read 
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For our correlated fits we have determined the covariance matrix according to the rules 
provided in statistics review of Ref. [13] treating the 10% uncertainty in the contribution of 
the continuum region for y/s > 11.21 GeV as a single independent error source. Explicit 
formulae for the covariance matrix can also be found in Ref. [14j . 
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